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ABSTRACT 

We study the properties of Lyman-a emitters (LAEs) and Lyman- break galaxies (LBGs) at z = 3 — 6 
using cosmological SPH simulations. We investigate two simple scenarios for explaining the observed 
Lya and rest-frame UV luminosity functions (LFs) of LAEs: (i) the "escape fraction" scenario, in 
which the effective escape fraction (including the IGM attenuation) of Lya photons is fhya ~ 0.1 
(0.15) at z = 3 (6), and (ii) the "duty cycle" scenario, in which the fraction of LAEs that are turned 
on at z = 3 (6) is Cduty ~ 0.07 (0.2) after correcting for the IGM attenuation. Our comparisons with 
a number of different observations suggest that the duty cycle scenario is preferred over the escape 
fraction scenario. We find that the mean values of stellar mass, metallicity and black hole mass 
hosted by LAEs are all smaller in the duty cycle scenario than in the escape fraction scenario. In 
our simulations, the galaxy stellar mass function evolves rapidly, as expected in hierarchical structure 
formation. However, its evolution is largely compensated by a beginning decline in the specific star 
formation rate, resulting in little evolution of the rest-frame UV LF from z = 6 to 3. The rest- frame 
UV LF of both LAEs and LBGs at z = 3 & 6 can be described well by the duty cycle scenario provided 
the extinction is moderate, E(B — V") « 0.15, for both populations, although our simulation might be 
overpredicting the number of bright LBGs at z = 6. We show that the correlation function of LAEs in 
the duty cycle scenario has a smaller correlation length than in the escape fraction scenario. The bias 
of LAEs relative to the dark matter distribution is higher at z = 6 than at z = 3 in both scenarios. 
The Lya LFs at z = 6 in a field-of-view of 0.2 deg 2 show a significantly larger scatter owing to cosmic 
variance relative to that in a 1 deg 2 field, and the scatter seen in the current observational estimates 
of the Lya LF can be accounted for by cosmic variance. 

Subject headings: cosmology: theory — stars: formation — galaxies: high-redshift — galaxies: evolu- 
tion - galaxies: formation - methods: numerical 



1. INTRODUCTION 

Hydrogen is ubiquitous in our universe, and its Lya 
emission line is now co mmonly observed in the spectra 
of high-z galaxies (e.g.. IShaplev et "al1l2003t lOuchi et al.1 
I2005U IHu fc Cowiel 120061). Recent detections of LAEs 
in large numbers at z = 3 — 6 using narrow-band fil- 
ters have opened up a new window for studying high-z 
galaxies after the long effo rt since the original proposal 
([Partridge fe Peebleslll967l ) of using Lya emission for the 
search of high-z galaxies. 

In particular, wide-field surveys have been ex- 
tremely successful in searching for LAEs, and the 
sizes of observed samples are becoming comparable 
to that of L BGs w it h several h undre d sources (e.g., 
Ouchi et all l2003afc fHuet alJ 12004 iDawson et all 
2004 iMalhotra fe Rhoadsl l2004t lOuchi et all 
2005al iTaniguchi et alJ 120051: IVenemans et a lJ 1 20051: 
Shimasaku et al ! I2006h __ Kashika waet all [2006; 
Gronwall etTaD 120071 : iMuravama et all 120071 ). These 
large samples allow us to construct the luminosity 
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function (LF) of LAEs with reasonable accuracy and to 
perform statistical analyses and comparisons with other 
populations such as LBGs, DRGs, EROs, BzKs, etc. 

At the same time, deep observations of smaller samples 
at infrared (IR) wavelengths are beginning to constrain 
the physical properties of LAEs such as st ellar mass and 
star fo rmation rate (SFR). For example, iGawiser et al.l 
(|2006h stacked the SEDs of 18 LAEs at z ~ 3.1, and es- 
timated an average stellar mass of M* ~ 5 x 10 s M Q 
and a star formation rate of SFR ~ 6 ft.f 2 M yr _1 
by fitting the SEP u sing a population synthesis model. 
IGawiser et al.1 (|2007f) updated the result using 162 LAEs 
from the MUSYC survey and reported M* ~ 1.0±°a x 
10 9 M Q and SFR ^ 2 ± l M^yr" 1 for typical LAEs at 
z = 3.1. lLai et all (|2007bl ) analyzed the same sample of 
z = 3.1 LAEs supplemented by the Spitzer IRAC ob- 
servations, and found M* = 3j^j x 1O 8 M , an age of 
~ 200 Myr, and an average SFR of 2M yr~ 1 for the 
IRAC(3.6^m)-undetected sample. The IRAC-detected 
sample is more massive with M* = 9 ± 3 x 1O 9 M , 
age = 1.6 ± 0.4 Gyr, an d SFR ~ 6 Mp, yr " 1 . 

At higher redshifts, iPentericci et alJ (|2007t ) studied 
z ~ 4 LBGs with and without Lya emission, and con- 
cluded that the LBGs with Lya emission are on average 
much younger and less massive than the LBGs without 
Lya emission. They estimated M* = (5 ± 1) x 1O 9 M 
and a n age of 200 ± 50 My r for the LBGs with Lya emis- 
sion. iPirzkal et all (|2007l ) studied the SEDs of 9 LAEs 
in the Hubble Ultra Deep Field at 4.0 < z < 5.7, and 
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estimated = 10 6 - 10 8 M Q and SFR w SA^yr" 1 . 
lLai et all ()2007af ) employed near-IR data of the Spitzer 
IRAC to increase the reliability of stellar mass estimates, 
and derived M* = 10 9 - 10 10 M© and ages of 5 - 100 Myr 
using 3 LAEs at z ~ 5.7 in the GOODS northern field. 
Overall, these observations suggest that LAEs are in gen- 
eral less massive, have lower SFRs, and are younger than 
LBGs. 

It has been suggested that LAEs might be ex- 
hibiting a very early phase of galaxy formatio n (e.g., 
IHu fc McMahonl 119961 : iMori fc Umemural I2006D . where 
the Lya photons emitted from the photoionized gas 
around star-forming regions are still able to escape from 
a relatively dust-free environment. However, it is also 
expected that Lya emission m ight be obseryed fr om a 
later stage of galaxy formation (jShaplev et al .1120031 ) . be- 
cause the galactic outflows driven by SNe can evacuate 
the dusty gas around star-f orming regions, as observed 
in loca l starburst galaxies fH eckman 2 0011) and high- z 
LBGs (jPettini et afl[2^ . [200l I Adelberger et~alll2003] ) . 

These ideas suggest evolutionary transitions from an 
early LAE to a LBG, and then back to a LAE at a 
later time. If this scenario is correct, then LAEs could 
be short-lived phenomena that occur only at certain 
phases of galaxy formation, and we need to consider 
the duty cycle of LAEs, similar to that of quasars (e.g., 
lHaiman fc H ul l200l. 

Furthermore. iKashikawa et al.l (|2007t ) detected an in- 
teresting segregation between LBGs and LAEs around a 
high-z quasar, and proposed that the existence of LAEs 
might be suppressed by the intense UV radiation from 
the nearby quasar. This suggests that the existence of 
LAEs could be affected by various parameters of the lo- 
cal environment, such as overdensity or intensity of the 
UV radiation field. 

Another interesting development in recent observations 
is that the LF in the rest-frame UV continuum of LAEs 
begins to be constrained at the same time, as well as the 
Lya LF. This sets important additional constraints to the 
evolution of LAEs from z = 6 to 3, requiring the models 
to fit the UV LF of LAEs as well as the Lya LF. The UV 
LFs of LBGs at z = 3 — 6 are being measured with in- 
crcasmg accuracy (e.g., Ouchi et al. 2004a; Yoshida et al. 
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20071) . and the comparison between the UV LFs of LBGs 
and LAEs would give us information about the relation- 
ship between LBGs and LAEs, which is one of the key 
questions that is intensely studied in the field. 

Motivated by these observations, there have been sev- 
eral theoretical studies of LAEs using semianalyti c mod- 
els of galaxy formation. iLe Delliou et alJ (120051 l2006h 
used the semianalytic model of iBaugh et al.l (|2005( ) to 
show that the Lya LF of LAEs at z — 3 can be explained 
if a uniform escape fraction of Lya photon / esc = 0.02 
is assumed. iKobavashi et al.l (120071) used the semiana- 
lytic model of iNagashima et all (|2005f ) and introduced a 
variable escape fraction of Lya photons owing to galactic 
wind feedback in high-z galaxies. 

In this paper, we use cosmological hydrodynamic sim- 
ulations of galaxy formation based on the concordance 
A cold dark matter (CDM) model to study the prop- 
erties of LAEs, with a focus on the evolution of SFR, 
Lya fc UV LFs, stellar mass and clusterin g. Pre viously, 
iNagamine et al.1 (|2004cft and iNight et all (|2006l ) exam- 



Run 


^box 


N p 


moM 


^gas 


€ 


Q5 


10. 


2 X 324 3 


2.12 x 10 6 


3.26 x 10 5 


1.23 


Q6 


10. 


2 x 486 3 


6.29 x 10 5 


9.67 x 10 4 


0.82 


D5 


33.75 


2 x 324 3 


8.15 x 10 7 


1.26 x 10 7 


4.17 


G6 


100.0 


2 x 486 3 


6.29 x 10 8 


9.67 x 10 7 


5.33 



Note. — Simulations employed in this study. The initial 
number of gas particles is equal to that of dark matter par- 
ticles, hence x2 for Np. the masses of dark matter and gas 
particles (moM and m gas ) are given in units of Ii~ 1 Mq, 
respectively, and e is the comoving gravitational softening 
length in units of kpc. The value of e is a measure of 
spatial resolution. All runs adopt our "strong" wind feed- 
back model. 

ined the properties of LBGs at z = 3 — 6 using the 
same set of cosmological simulations analyzed here, and 
showed that the rest-frame UV LF of LBGs can be ex- 
plained reasonably well with moderate extinction. Build- 
ing on top of our previous work on LBGs, here we discuss 
the relationship between LBGs and LAEs. 

We have so far only discussed the point sources 
of Lya emission, but there are also source s with ex- 
tended Lya emission, called "Lya blobs" dKeel et al.l 
[19991 ISteidel et al.l l2000t IMatsuda et"aLll2004D . One of 
the possibilities is that the extended Lya emission is 
powered by the release of gravitational potential en- 
ergy as the baryons condense inside d ark matter halos 
(jHaiman et al.ll2000l : iFardal et al.ll200ll ) Using the same 
series of simulations as in this paper, iFurlanetto et all 
(2005) considered the Lya emission from diffuse IGM 
and the gas in halos around galaxies, which is powered 
by gravitational processes and photoionizing background 
radiation. They concluded that the Lya emission from 
recombinations that follow the absorption of stellar ion- 
izing photons (i.e., the Lya emission associated with star 
formation) dominate s the total L ya ph oton production 
rate. Furthermore, iSaito et al.l (|2007D found that the 
number density of Lya blobs are only 10 — 20% of the 
total LAE population. Therefore, in this paper, we focus 
on the Lya emission associated with star formation and 
do not consider Lya blobs. We concentrate on the com- 
parison with updated observations of LAE LF and their 
relationship to those of LBGs. 

This paper is organized as follows. In §[^1 we briefly de- 
scribe our simulations. We then discuss the specific SFR 
and stellar masses of high-z galaxies in our simulations in 
§ [3] We present the evolution of the galaxy stellar mass 
function in § [H and discuss the evolution of the Lya LF 
from z = 6 to 3 in § [5] We propose two simple scenarios 
to explain the Lya LF of LAEs. For both scenarios, we 
compute the mean values of stellar mass (§0, black hole 
masses (§ and metallicity (§ [5]) of LAEs. We discuss 
the relationship between Lya and UV LFs in § [51 and 
their evolution in § [TDJ The correlation function of LAEs 
is presented in § 1111 Finally, we discuss the cosmic vari- 
ance in the current surveys of LAEs in § fTS] and conclude 



2. SIMULATIONS 



We use the smoot hed particle hydrodynamics 
(SPH) code GADGET2 (|Springell I2005D in this work. 
It employs the 'entropy conserving' formulation 
(Springel fc Hernquist 2002) to alleviate the overcooling 
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Fig. 1. — Specific star formation rate vs. galaxy stellar mass at z = 3 and 6. The three sets of contours are for the Q5 (black), D5 (blue) 
and G6 (red) run from left to right. The median in each mass bin is shown with the symbols. The analytic fits to the median points are 
shown with the long-dashed line and are given in the text. Galaxies with zero SFR are indicated at logSFR/M* = —4. 



problem, which previous generations of SPH codes 
experienced. Our simulations include radiative cooling 
by hydrogen and helium, heating by a uniform UV 
background (e.g., iKatz et alj 119961 : iDave et all I1999D . 
star formation and supernova fee dback based on a sub- 
partic le multiphase ISM model (|Springel k. Hernauistl 
l2003al ) , and a phenomenological model for galactic 
winds (|Springel fe Hernauistl l2003bl) . 

The details of the star for mation model were described 
in iNagamine et al.l (|2004bf) , so we only give a brief de- 
scription here. In short, gas particles are allowed to 
spawn a new star particle when a set of criteria is satisfied 
at each time-step. Groups of star particles are regarded 
as galaxies in the simulation, and we identify them by 
apply a grouping algorithm. The code tags the star par- 
ticles with physical quantities such as their mass, forma- 
tion time, and metallicity. Using these tags, we compute 
the spectrum of each star particle with th e population 
synthesis code of iBruzual fe CharlotJ (|2003D . and co-add 
the individual luminosities to obtain the spectra of our 
simulated galaxies. 

We use four different simulations with varying box sizes 
and particle numbers (see Table 1) in order to cover a 
wide range of halo masses and assess the resolution ef- 
fect. These simulations exte n d the set of runs carried 
out by iSpringel fe Hernauistl (|2003bD to higher resolu- 
tion. Unfortunately, the Q6 simulation was stopped at 
z ~ 4 owing to its very long computing time. Therefore 
we basically use the Q5 run in this paper, and show the 
results from the Q6 run for z = 6 where appropriate. The 
results of the Q5 and Q6 runs are very similar at z = 6, 
except that the Q6 run has a slightly better coverage 
for the lowest mass galaxies with M* < 1O 7 ' 5 M . The 
main conclusions of this paper are not affected by the ab- 
sence of Q6 results at lower redshifts. The adopted cos- 
mological parameters of all simulations considered here 
are (fi m , fi A , fib, 08, h) = (0.3,0.7,0.04,0.9,0.7), where 
h = ffoAiOOkms^Mpc- 1 ). 



3. SPECIFIC SFR AND GALAXY STELLAR MASS 

We start by examining the instantaneous SFR of sim- 
ulated galaxies at z = 3 and z = 6, because we will use 
this quantity to calculate the Lya luminosity. Figure [1] 
shows the specific SFR (= SFR/M*; SFR per unit stel- 
lar mass) as a function of galaxy stellar mass at z = 3 
and 6. The three sets of contours are for the Q5 (black), 
D5 (blue) and G6 (red) runs, from left to right. The 
Q6 result at z = 6 is very similar to that of Q5, so it is 
not shown here. Each simulation box can resolve only 
a limited range of galaxy masses, so we use three differ- 
ent simulations to cover a wide range of stellar masses, 
ill* = 10 7 - 1O 12 M . The median value of the specific 
SFR is shown by the symbols for each bin of log M*. The 
number of galaxies with no star formation is greater at 
z = 3 than at z = 6, as indicated at logSFR/M* = —4.0. 

The distribution broadens at the lower mass end of 
each contour for two reasons. One is that there is a 
larger number of lower mass halos, therefore the distribu- 
tion naturally becomes broader as the larger population 
exhibits a larger variation in its properties. The other 
reason is that the resolution limit of each run progres- 
sively shifts to lower masses, and close to the resolution 
limit the distribution broadens owing to numerical noise. 

At z = 3, the specific SFR is almost constant (with sig- 
nificant scatter around the mean) across the mass range 
of M* = 10 8 - 10 M Q . On the other hand, at z = 6, the 
specific SFR is an increasing function of galaxy stellar 
mass, indicating that star formation is more efficient in 
more massive galaxies. Star formation becomes rapidly 
inefficient in low mass galaxies with M* < 10 8 M Q , and 
the distribution seems to drop off completely at M+ ~ 
1O 7 M . This rapid fall-off of the SFR at M* ~ 10 7 M Q 
may be related to the threshold density for star formation 
in the simulati on and the obser ved SF cutoff in nearby 
spiral galaxies (jKennicuttl IT998). A detailed investiga- 
tion of the SF threshold in the simulation is beyond the 
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Fig. 2. — SFR function at z = 3 (thin lines) and z = 6 (thick 
lines) for the Q5 (long dashed), D5 (short dashed), and G6 (solid) 
runs. 



scope of this paper, and is deferred to future work. 

We provide the following approximate fit to the median 
points shown in Fig. [TJ 



Y = aX - 



X 



(1) 



where (X, Y) = (log M* , log SFR/ M* ) , and (a,b,c,d) = 
(0.0, 0.17, 7.0, -0.10) & (0.15, 0.2, 7.0, -0.75) for z = 3 & 
6, respectively. At z = 3, the down-turn at logM* « 7.0 
is not clearly seen, but we kept the value of c the same 
for both redshifts for simplicity. 

Figure [2] shows the galaxy SFR function at z — 3 and 6, 
which measures the differential number density of galax- 
ies per logarithmic bin of SFR, similarly to a galaxy lu- 
minosity function. We find that there is not much evo- 
lution in the SFR function from z = 6 to 3. This might 
seem counter-intuitive given the evolution seen in the 
specific SFR, however, the evolution of the SFR function 
is caused by a combined evolution in both the specific 
SFR and galaxy stellar mass functions. 

4. GALAXY STELLAR MASS FUNCTION 

In our simulations, the galaxy stellar mass function 
evolves rapidly from z = 6 to z — 3, as expected in a 
hierarchical universe (Figure [3]); new halos of low mass 
are constantly formed, less massive galaxies merge into 
more massive systems, and the number density of mas- 
sive galaxies increases with decreasing redshift. We note 
that galaxies grow in stellar mass in our simulations, but 
the star formation becomes less efficient from z — 6 to 3 
(Fig. [I]). The latter effect compensates for the growth of 
the mass function, resulting in little change in the SFR 
function (Fig.©. 

At z = 3, the s imula tion agrees well with the data 
from iDrorv et all (120051 . blue circles, 3.0 < z < 4.0), 
Fontana et al.1 (120061 green triangles, 2.5 < z < 3.0) and 
Perez- Gonzalez et al.l ( 20071 red squares, 3.0 < z < 3.5) 
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Fig. 3. — Galaxy stellar mass function at z = 3 (thin lines) 
and 2 = 6 (thick lines) for Q5 (long dashed; Q6 for 2 = 6), D5 
(short dashed), and G6 (solid lines) runs. The dotted line indi- 
cates the power-law of n(M*) oc M^" 2 2 . The data po ints are from 
IDrorv et aLl J2005I. blue circles) , IFontana et alj pOOti green trian- 
gles) and [Perez-Gonzalez ct al. (2007| red squares). 



at log M* > 10.5. The data by IDrorv et~aT1 (|2005l ) 
suggest that the mass function becomes shallower at 
logAf* < 10.0, but our simulations have a steeper slope 
(ra(M*) oc Afr 2 ' 2 ) at the low -mass end as indicated by 
the dotted line. While it appears likely that our simu- 
lations overpredict the number of low-mass galaxies, fu- 
ture deeper observations are needed to check this, based 
on a more reliable measurement of the faint-end of the 
mass function. The location of the "knee" in the simu- 
lated mass function is uncertain, given the steep faint-end 
slope. 

We compute the total stellar mass density by integrat- 
ing the interpolated mass function in Fig.[3]over the mass 
range of logM* « [7.0, 12.0], and obtain = 3.5 x 10 s 
(8.8 x 10 7 ) M Q Mpc" 3 at z = 3 (6). This corresponds 
to fl+ = 0.0026 (6.5 x 10~ 4 ) at z = 3 (6). These values 
ar e higher than any of the cosmic SFR models presented 
bv lNagamine et all (|2006bL Fig. 6a), suggesting that the 
faint-end slope in Fig. [3] might be too steep. Neverthe- 
less, we will use these total stellar mass densities in § [6] 
to compute the fraction of stellar mass that LAEs con- 
tribute. 

5. Lya LUMINOSITY FUNCTION 

We compute the intrinsic Lya luminosity L^ a emitted 
by high-z galaxies as 

V&a = W 42 (SFR/M yr- 1 ) ergs" 1 , (2) 

following iFurlanetto et al.l (|2005h . This relationship is 
accurate to within a factor of a fe w according to the stel - 
lar population synthesis model of iLeitherer et all (1999) 
for a Salpeter initial mass function (IMF) with a mass 
range of 1 — 100 M Q and metallicities between 0.05 < 
Z/Zq < 2. We use the simulated SFR in the right-hand- 
side of Eq. (0). 
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Without any corrections to the Lya luminosity, we find 
that our simulations overpredict the Lya LF by a sig- 
nificant factor (~ 10) compared to the observations of 
lOuchi et all IpOOl . H ere, we choo s e to c ompare our re- 
sults with the data bv lOuchi et al.l (|2007f h because their 
sample comes from a large survey field and they have 
performed extensive comparisons with the earlier LF es- 
ti mates. Ouchi et al.'s L F at z ~ 3 is consistent with that 
of iGronwall et all (|2007| ). According to their data, there 
is not much evolution (no more than a factor of 2 — 3) 
between z — 6 and 3 in the observed apparent Lya LF 
either in luminosity or number density. However, Lya 
fluxes from high-z sources are attenuated by the inter- 
galactic neutral hydrogen, causing an asymmetric profile 
in the Lya emiss ion line with the blue-side being ab- 
sorbed more (e.g..lHu et aI1l2004t iKashikawa et a,l.ll2006t 
IShimasaku et al.ll2006h . Therefore, when the data is cor- 
rected for this effect, little evolution in the apparent Lya 
LF means strong evolution in the intrinsic LF, in the 
sense that the Lya luminosity and/or the number den- 
sity of LAEs are intrinsically brighter/higher at z ~ 6 
than at z ~ 3. 

In the following, we consider two possible scenarios to 
match the simulation results to the observed apparent 
Lya LF. The two proposed scenarios are very simple, but 
they capture the two extreme situations that plausibly 
bracket the true behavior. 

5.1. Escape Fraction Scenario 

In the first scenario we simply assume that only a fixed 
fraction of Lya photons reaches us from the source, i.e., 

pobs r ^intrinsic /o\ 

Lya ~ iLytt^Lya ! \°) 

where Fh ya is the Lya flux. The parameter fhya can 
be interpreted as an effective escape fraction that in- 
cludes the following three effects: escape of ionizing pho- 
tons, local dust ext inction, and absorption by the IGM 
([Barton et al.ll200^ ). We characterize this as 

fhya = /dust (1 ~ / C sc) /iGM) (4) 

where /dust is the fraction of Lya photons that is not 
extinguished by local dust, / ™ is the fraction of ionizing 
photons that escape from galaxies and thus create no Lya 
photons, and /igm is the fraction of Lya photons that 
are not absorbed by the IGM, i.e., the transmitted flux. 
We call this case the "escape fraction" scenario. 

Of course, in the real universe, different galaxies may 
have different values of /dust and f%£?, depending on their 
physical parameters such as age, mass, SFR and local 
environment. These parameters can also depend on red- 
shift. Therefore, the above parameterization should be 
interpreted as an attempt to capture the average behav- 
ior of bright galaxies that are currently being observed, 
even though for simplicity we do not indicate the implicit 
averaging with (• • ■} in our notation. 

The left column of Figure [4] shows a comparison of 
our simulat i on re sults with the observational data by 
lOuchi et alJ (|2007l ). adopting / LyQ = 0.1 (0.15) for z = 3 
(6). This scenario corresponds to simply shifting the 
simulated LF toward lower luminosity, therefore the cur- 
rently observed LAEs correspond to relatively massive 
galaxies with high SFR. Here we selected the values of 
fhya such that the G6 run agrees well with the observed 



data points, because this run has the largest box size 
and covers the bright-end of the observed LF much bet- 
ter than our other runs. The D5 run underestimates the 
number density of massive galaxies with logLLya ^ 42 
owing to its smaller box size. The agreement between the 
simulation results and the observed data is very good at 
both z = 3 and 6, including the slope of the LF. Since 
our SFR function does not evolve very much (Fig. [2]), the 
values of fhya at z — 3 and 6 are very close. 

5.1.1. 5.1.1. IGM attenuation and /J™ , /dust 

We estimate the effect of IGM attenuation to be 

/igm = e- Trff = 0.82 (0.52) for z = 3 (6) (5) 

using the iMadaul (1995) formulation with the assump- 
tion that only half of the symmetric Lya line is ab- 
so rbed. These values are consistent with those obtained 
bv lOuchi etaH (|2007f ). Inserting Eq. © into Eq. (gj), we 
obtain 

/dust (1 - Cc) = 0-12 (0.29) for z = 3 (6). (6) 

IChen et ahl (|2007f ) reported that, using the afterglow 
spectra of long-duration gamma-ray bursts, the mean es- 
cape fraction of ionizing radiation from sub-L* galaxies 
at z > 2 is (Z™ 1 ) = 0.02 ± 0.02 with an upper limit of 
(/esc ) — 0-075. If the escape fraction of ionizing pho- 
tons is as small as /g°" = 0.02, then our result implies 
i riust ~ 0-12 (0 . 29) at z = 3 (6). 

Ilnoue et alJ (|2006l Fig. 3) compiled existing direct 
measurements of escape fractions of ionizing photons and 
estimates based on the observed ionizing background in- 
tensities. They suggested that the value of /™ might 
be increasing with redshift: /j°° w 0.02, 0.06 & 0.2 at 
z = 2, 3 & 4 — 6. In this case, 

/ dust = 0.13 (0.36) at z = 3 (6). (7) 

The lower value of /dust at z = 3 suggests that the envi- 
ronment around the star-forming regions becomes more 
polluted by dust as star-formation proceeds from z = 6 
to 3, blocking more Lya photons. 

Many researchers (e.g., iLe Delliou et al.1 120061 : 
iKobavashi et alJ l2007h simply adopted /igm = 1-0, 
arguing that various effects can reduce the amount of 
IGM attenuation, such as ionization of the IGM around 
galaxies, clearing of the IGM by galactic winds, and 
rcdshifting of Lya photons by the scattering in the 
wind. They also refer to the fact that the reionization 
of the Universe was mostly completed by z ~ 6, as 
indicated by measurements of Gunn-Peterson absorption 
in quasar spectra, and a s suggested by constr aints on 
the clustering of LAEs (|McQuinn et all [2007] ). It is 
possible that the asymmetric Lya line profile is caused 
by the local ISM at the source, rather than by the IGM. 
Given the large uncertainty in the value of /igm, we 
also consider the case of /igm = 1-0. In this case, our 
result implies / dust (1 - = 0.10 (0.15) for z = 3 (6). 
Ado pting the values of f™ = 0.06 (0.20) at z = 3 (6) 
from Ilnoue et al.l (|2006h . we obtain 

/ dust = 0.11 (0.19) at z = 3 (6). (8) 

In either case, our results imply / du st ~ 0.1 at z = 3, 
and /dust ~ 0.2 — 0.4 at z = 6. 
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5.2. Duty Cycle Scenario 

The other scenario we examine is based on the assump- 
tion that only a fixed fraction (Cduty) of all galaxies is 
active as LAEs at a given time. As far as the luminosity 
function is concerned, this is equivalent to assuming that 
each LAE has a certain duty cycle and remains turned 
on only for a fixed duration of time, therefore, we call 
this case the "duty cycle" scenario. 

Before we calculate the values of Cduty, we correct the 
observed Lya LF for the effect of IGM attenuation by 
factors of (/igm)^ 1 = (0.82)" 1 & (0.52)- 1 for z = 3 & 6, 
respectively. As shown in the right column of Figure 01 



we obtain good agreement with the data if we assume 
Cduty ~ 0.07 (0.2) for z — 3 (6). This is equivalent to 
lowering the normalization of the simulated LF to match 
the observed data. Here we adjust our simulated LF so 
that the results of the D5 and G6 runs bracket the ob- 
served data points, because in this scenario the observed 
data overlaps with lower mass galaxies in the D5 run. 
The above value of Cduty can be interpreted as either 
only 7% of the sources are turned on as LAEs at z = 3, 
or LAEs are turned on only for 70 Myrs out of 1 Gyr at 
z«3. 

The yellow shaded region in the lower panels of Fig- 
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ur e [H indicates t he re gion covered by the data points 
of iSantos et al . (2004) and the simulation results of 
iDave et al.1 (|2006T) . Their results are significantly lower 
than our simulation results at log Llvq = [40.5, 42.5] . We 
comment further on this point in §[l2] and § 1131 

6. STELLAR MASS OF LAES 

Figure [5] summarizes the relationship between galaxy 
stellar mass, SFR and Lya luminosity. The yel- 
low shaded regions indicate the luminosity range of 
logLLya = [42,44], which roughly corresponds to th e 
currently observed LAEs at z = 3-6 (jOuchi et al.ll2007l ). 
The corresponding stellar mass ranges differ significantly, 



depending on the two scenarios and redshift as summa- 
rized in Table [2j but they are not affected by the un- 
certainties in the value of /igm, because Fig. [5] is solely 
determined by Eq. ([2]) and the values of fhya- 

In the escape fraction scenario, the raw (i.e., before 
any corrections) simulated Lya LF is simply shifted to- 
ward lower luminosity without a change in the normal- 
ization, therefore currently observed LAEs correspond to 
the most massive and luminous objects at the brightest 
end of the LF. The mean stellar masses of LAEs with 
log-Lxya = [42,44] in the G6 run are 

(M^) = 2.5 x 10 10 (1.9 x 10 9 ) M Q at z = 3 (6), (9) 
respectively, as given in Table [2j These values are close 
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Fig. 6. — Galaxy stellar metallicity vs. galaxy stellar mass. The three sets of contours are for the Q5 (black), D5 (blue) and G6 (red) 
run from left to right. For z = 6, the result of the Q6 run is also shown with the green contour. The range of galaxy stellar mass for each 
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to the lower edge of the listed stellar mass ranges, be- 
cause of the increasing number of low-mass galaxies in 
a cold dark matter universe. In this scenario, the LAEs 
with logi LyQ = [42,44] contribute only 18% (11%) of 
the total stellar mass density at z = 3 (6). This rela- 
tively low fraction owes to the steep faint-end slope of 
the stellar mass function (Fig. [3]). 

In the duty cycle scenario, currently observed LAEs 
correspond to slightly lower mass galaxies than in the 
escape fraction scenario: 

(M*) = 3.8 x 10 9 (6.1 x 10 8 ) Mq at z = 3 (6) (10) 

for the LAEs with logL Lya = [42, 44] in the D5 run. The 
fraction of stellar mass density contributed by the LAEs 
is higher in this scenario, amounting to 42% (29%) for 
z = 3 (6). 

7. BLACK HOLES HOSTED BY LAES 

The masses of black holes hosted by LAEs are also 
of significant interest, as they determine the AGN con- 
tribution to the total energy output from LAEs. We 
indicate the black hole masses in the top axes of Fig- 
ureEJassuming M*/M BH ~ 0.004 (0.005) at z = 3 (6), as 
suggested by the recent numerical simulations of galaxy 
mergers jD i Matt ecTet al.l [2005; Ro bertson et alJ 120061: 



Hopkins et al.ll2007ch and observations (e.g. iPeng et ail 
2006J). We find that the LAEs with logL Ly a = [42,44] 



host BHs with masses log -Mbh = [7.7,9.7] (z — 3) and 
[6.9,8.7] (z — 6) for the escape fraction scenario, and 
log Mbh = [6.7,8.8] (z = 3) and [6.3,8.0] (z = 6) for 
the duty cycle scenario. If we instead assume the lo- 
cal v alue of M*/M B h ~ 0.001 (jKormendv fc Gebhardtl 
1200 ll ) , then the above mass ranges change to log Mbh = 
[7.2,9.2] (z = 3) and [6.3,8.1] (z = 6) for the escape 
fraction scenario, and logA/eH = [6.2,8.3] (z — 3) and 
[5.7, 7.4] (z = 6) for the duty cy cle scenario. 

At z = 3.1. IOuchi et alJ (|2007ft reported that the AGN 
fraction of LAEs is only ~ 1% for LAEs with L^ ya > 



1 x 10 42 ergs -1 , but that the fraction is as high as 100% 
for the bright LAEs with L Lya > 4 x 10 43 ergs -1 . On the 
other hand, they found no AGNs (nor LAEs) at z = 5.7 
with a bright luminosity of L^ ya > 4 x 10 ergs -1 , 
and suggested that the number density of LAEs with 
AGN activities would drop from z = 3 to 6. Fig- 
ure [5] shows that the mean BH mass for the LAEs with 
log LLya = [42,44] decreases from z = 3 to 6, which is 
at least consistent with the observed decrease of bright 
AGNs if BH mass is correlated with the strength of AGN 
activity. 

8. LAE METALLICITY 

In Figure El we show the galaxy stellar metallicity vs. 
stellar mass in our simulations. The metallicity of each 
galaxy is computed by summing up the metal mass in all 
constituent star particles, and then dividing by the total 
stellar mass of each galaxy. There is a weak positive 
correlation between the two quantities with significant 
scatter. 

At z = 6, the majority of galaxies have \og(Z/Z Q ) < 
—0.5, and the most massive ones with logM* > 10.5 
have \og(Z / Zq) > —0.5. There are some outliers with 
log(Z/Z Q ) > —0.5 at the low-mass end of the distribu- 
tion, which are small galaxies that have just started to 
undergo star formation and self-enrichment. Close to the 
resolution limit, the coarse sampling of the galactic winds 
leads to substantial numerical scatter in the metal loss, 
which in turn can temporarily produce high metallicities 
for some of the small galaxies. However, the fraction of 
such outliers is very small with ~ 6% (0.2%) of the total 
sample in the Q5 run at z — 3 (6). 

At z = 3, the majority of galaxies still has log(Z/Z©) < 
—0.5, but the number of galaxies with log (Z/Z©) > — 0.5 
has increased significantly since z — 6, and the most mas- 
sive ones with logM* > 11 approach solar metallicity. 

Since the galaxy stellar mass range is lower in the duty 
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cycle scenario than in the escape fraction scenario, the 
corresponding galaxies have lower metallicity in general, 
as summarized in Table [2] Owing to the large scatter in 
the distribution, the difference in the metallicity range 
is not so large between the two scenarios, but the trend 
of lower metallicity in the duty cycle scenario is clearly 
seen in the mean metallicity values listed in Tabled As 
expected, the mean metallicity at z = 3 is higher than 
at z = 6 by about a factor of 2 in both scenarios. 

It would be useful to compare the mean metallicity 
of LAEs and LBGs. The mean metallic ity of LBGs is 
known to be ~ 1/3 Z e (e.g.. lPettmll2004l ). At z = 3, the 
mean metallicity of LAEs in the escape fraction scenario 
is comparable to that of LBGs with (Z/Zq) = 0.39, while 
it is lower for the duty cycle scenario with (Z/Zq) = 0.21. 

9. REST-FRAME UV LUMINOSITY FUNCTION 

One of the interesting recent development in the obser- 
vations of LAEs is that the rest-frame UV LF of LAEs is 
beginning to be constrained at the same time as the Lya 
LF. This has been difficult in the past, because the UV 
continuum of LAEs tends to be faint, and large samples 
of LAEs were not available owing to limited sizes of the 
field-of-view (FoV) of the observations. 

Figure [7| compares the simulated rest-frame UV LF 
with observational data. We adopt a uniform, moder- 
ate extinction of E(B — V) = 0.15 at both z = 3 and 
6 following our p revious work (Nagaminc "eT^l2004cl 
iNight et al.ll200a K This exti nction value is the median 
value for the LBGs at z ~ 3 (jShaplev et al.ll2001h . 

At z — 3, we obtain very good agreement between 
the G6 run (blue sol i d line ) and the observed data of 
lAdelberger fc Steidell (2000, open squares) for LBGs. 
This agreement suggests that the SFR function (Fig. [2]) 



of simulated galaxies at z = 3 is quite reasonable at least 
at the bright-end. We note that it is also reasonable that 
the D5 run underpredicts the observed data owing to its 
limited box size. 

At z = 6, we show in Fig.[7]the data bv lBouwens et al.l 
(|2007| ) for z-dropout LBGs. The yellow shade indicates 
the region covered by the Schechter function fits of other 
obser v ational studies (|Bpuwens et al.ll 2004: Bunke r~et al.l 
2004 IDickinson et al l l2004t lYan fc Windhorst! 12004 
Malhotra et al.l [200 5; Beck with et alJ 120061), wh ich we 
estimated from Fig. 11 of iBouwens et al.l (|2007D . The 
G6 run follows the upper envelope of the yellow shade 
and slightly q v erpred icts the data at Mjjv < —22. 
IBouwens et al.1 (|2006t ) reported that the extinction of 
z = 6 LBGs are lower than those of z = 3 ones, which 
would exacerbate the discrepancy between the simulation 
and observation at z — 6. 

It is possible to improve the agreement between the 
G6 run and the observed data by increasing the assumed 
extinction for massive galaxies. For example, correlating 
the extinction with metallicity would boost the extinc- 
tion in massive galaxies (Fig. [6]) , and make the bright-end 
of the LF steeper. However, here we choose not to resort 
to such an additional modeling, because this would cloud 
the interpretation of the comparison between simulation 
and observation. 

Another possible cause for the overprediction of UV LF 
at the bright-end at z = 6 is that the current simulations 
lack the explicit implementation of AGN feedback. It is 
considered that the energy and momentum feedback from 
black holes suppress the star form ation in massive galax- 
ies at intermediate redshifts (e.g..lDi Matteo et al.ll2005l : 
ISpringel et aLll2005h iHopkins et al.ll2007bl )ah . The effects 
of AGN feedback on galaxy LFs at high redshifts can 



10 



Nagamine et al. 



be studied in the future using hydrodyn amic simulations 
that treat AGN feedback explicitly (|Di Matteo et al.1 
[2007tlSiiacki et al.ll20"07t iLi et al.ll200l . 

10. EVOLUTION OF Ly« AND UV LF OF LAES 

Currently there is no clear picture regarding the evo- 
lution of the observed rest-frame UV LF from z = 6 
to 3. Some studies suggest a systematic brightening of 
the characteristic magnitude M^ v rest by 0.5 — 1 mag 
(IBouwens et al.1 120061: lYoshida et ail 120061 : lOesch et all 
l2007t IBouwens et al.ll2007j ). Oth er work argues for an 
evolution in the normaliz ation 4>* ( Bcck with et aLlfeOOGf) 
or m the faint-end slope (|Yan fc Windhorst) l2004f ). 

In our current simulations, there is not much evolu- 
tio n in the UV LF fr om z — 6 to 3, as already shown 
bv iNight et~aTl (|2006l ). except that the M T * TV becom es 
slightly brighter (see Fig. 13 of IBouwens et al.ll2007l in 
which the comparison between simulations and observed 
data was performed explicitly) and a becoming slightly 
shallower, as can be seen in Fig. [7]when compared to the 
power-law slope of a — —2.2 (dotted line). This trend is 
at least qualita t ively consistent with the observed one by 
IBouwens et al.l (|2007| ). The galaxy stellar mass function 
evolves as shown in Fig. [3j but the evolution in the SFR 
(Fig. (2) cancels that out, resulting in little evolution in 
the UV LF. 

Also shown in Fig. [7] are the UV LF of LAEs by 
lOuchi et al.l (|2007l filled triangles). Their data suggest 
that the number density of LAEs at z = 3 is only 10% 
that of LBGs', down to Muv — —20 for the equivalent 
width (EW) limit of 40 - 60 A, whereas the LAE frac- 
tion at z ~ 6 increases to 50 — 100% for the EW limit of 
~ 30 A. Ouchi et al.'s LAE fract ion at z = 3 i s con- 
sistent with the earlier result by IShaplev et alJ ()2003L 
~ 25% with EW limit of 25 A), but at z ~ 6, Ouchi 
et al.'s LAE fraction is higher than other spectroscopic 
studies of ^-dropout galaxies, which typically suggest 
~ 30% dDow-Hvgelimd et alJl2007t IStanway et a 1.1 1200 ll: 
IVanzella et al.ll2006fh 

Let us first consider the duty cycle scenario. At z = 3, 
when we multiply the normalization of the simulated 
UV LF (which assumes uniform E(B — V) = 0.15) by 
Cduty = 0.07 (i.e., LAE fraction of 7%), we obtain good 
agreement with the observed UV LF of LAEs. At z = 6, 
we again obtain a reasonable agreement with the ob- 
served UV LF of LAEs when we multiply the LBG UV 
LF by Cduty = 0.2 (i.e., LAE fraction of 20%), as shown 
in the right panel of Fig. [7] with the short-dash long- 
dashed line. However, given that the LF of z'-dropout 
galaxies at z ~ 6 is still very uncertain and that our 
simulated LF might be overpredicting the UV LF at the 
bright-end, we cannot make a strong argument about the 
LAE fraction at this point. If current observations are 
underestimating the number of very bright z'-dropouts 
(Myv S —22) at z — 6 for some reasons (e.g., cosmic 
variance, limited FoV), then our simulations would be 
able to explain both the UV and Lya LF of LAEs at 
z = 3 & 6 with LAE fractions of 7% & 20%, respectively, 
in the duty cycle scenario. 

One problem in the above argument is that we im- 
plicitly assumed that the UV extinction of LAEs is 
also E(B — V) = 0.15. However , rece nt observations 
(|Ouchi et alJ 120071 : iGronwall et~aT1 120071 ) suggest much 
lower extinction for LAEs with E(B — V) < 0.05. If this 



is the case, then the duty cycle scenario would not be able 
to explain the UV and Lya LF of LAEs self-consistently 
with above values of Cduty Furthermore, the escape frac- 
tion scenario completely fails to explain the UV LF of 
LAEs if the extinction of LAEs is systematically lower 
than that of LBGs. Because this scenario assumes that 
"all LBGs emit Lya emission following Eq. Q that are 
then attenuated uniformly by a factor of fhya" , the es- 
cape fraction scenario would predict a brighter UV LF 
for LAEs than that of LBGs, which is contradictory to 
the observed data. 

At this point there is no easy way out, unless we start 
considering the possibility that the LAEs are a separate 
population from LBGs. It is certainly possible that the 
majority of LAEs are low-mass galaxies with low extinc- 
tion and low metallicity. This is also consistent with what 
we have argued for in Fig. [5] for the duty cycle scenario. 
Even if we bring a fraction of low-mass galaxies in the 
UV LF with E(B - V) = 0.0 and multiply it by a cer- 
tain value of Cduty, the slope of the LAE UV LF does not 
agree with the data. This suggests variable extinction for 
LAEs, and we could model this by correlat ing extinction 
with metallicity (e.g.. iFinlator et "aT1l2006h . although we 
discuss the problems in this approach in § [T3] Since the 
purpose of our paper is not to introduce models more 
complicated than the simulation itself, we choose not to 
pursue any additional modeling. 

11. CORRELATION FUNCTION AND BIAS OF LAES 

The correlation function of LAEs provides interesting 
constraints on the distribution of LAEs, and it might 
help to discern the two scenarios that we have discussed 
above. We compute the auto-correla ti on fu nction (CF) 
of LAEs using the lLandv fc Szalavl (|1993f ) estimator, 
(DD-2DR+RR)/RR, and show the results in Figured 
Since the current luminosity limit of the observed LAEs 
is logLLya ~ 42.0 for the z — 3 sample, we restricted our 
sample to those with log-Lxya > 42.0 for both z = 3 & 6 
to measure the correlation function. At z = 6, the cur- 
rent luminosity limit is higher at log-Lxya > 42.5 because 
of the greater distances to the sources, but we kept the 
luminosity limit the same for our calculations to obtain 
a reasonably strong correlation signal with a sufficient 
number of simulated galaxies. 

In the escape fraction scenario (Fig. EH), the observed 
luminosity range is completely covered by the G6 run, 
so it is well-suited for this calculation. Using the data 
points in the range of logr = [0.0,1.3], we perform a 
least square fit to the power-law function £(r) = (r /r) 7 , 
and find {r [hT l Mpc], 7 ) = (5.5, 1.67) and (4.4, 1.68) for 
z — 3 & 6, respectively. The CF at z — 3 has a longer 
correlation length than at z = 6, because the mean stellar 
mass of LAEs at z = 3 is higher ( § EJ) - The slopes of the 
CF at the two epochs are very close to each other with 
7 w 1.7. Our val ues of (r n ,7) are similar to those of 
LBG s at z w 3 - 6 (|Ouchi et alj|2004ri lAdelberger et all 
I2005T) , which makes sense because the mean stellar mass 
of LAEs in the escape fraction scenario is comparable to 
that of LBGs with a few times 10 10 M Q at z = 3. 

For the duty cycle scenario, the observed data points 
of lOuchi et al.l (|2007l ) are fully covered only by the D5 
run (see Fig. |4|). The resolution limit of the G6 run 
corresponds to logLxya ~ 43 in this scenario, therefore 
we cannot use the G6 run to measure the CF for the 
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duty cycle scenario with the same luminosity limit of 
logLLya = 42.0. As discussed in § 15. 2\ the number den- 
sity of sources in the simulation has to be reduced by a 
factor of Cduty = 0.07 (0.2) at z = 3 (6) in this scenario. 
In the D5 run, there are 7V>42 = 1439 ( 20 3 2) sources 
with logLLya > 42.0 at z = 3 (6), therefore we need to 
select only iV>42 x Cduty = 101 (406) LAEs in the co- 
moving volume of (33.75 hr 1 Mpc) 3 . Owing to the small 
sample size, the CF signal for the duty cycle scenario 
is somewhat noisy. We therefore randomly resample 20 
different data sets with above LAE numbers, and cal- 
culate the mean of the 20 different trials to reduce the 
noise. We find that the CFs drop off at r > 8/i _1 Mpc 
owing to limited box size, therefore we only use data at 
r < 8 h^ 1 Mpc for the power-law fit. Fig. [8b shows that 
the CF at z = 3 is steeper with 7 = 2.3, although we con- 
sider that this result is not reliable owing to small sample 
size (101 LAEs). At z = 6, we obtain a shallower slope 
of 7 = 1.49 and ro = 3.1 /i" 1 Mpc. In fact, if we increase 
the sample size at z = 3 to the same size as that at z = 6 
(406 LAEs, corresponding to Cduty = 0.28, which will 
overpredict the Lya LF) , then we obtain a similar signal 
to that at z = 6 with 7 = 1.61 and ro = 3.9 h^ 1 Mpc. 
Therefore we consider that the steep slope of z — 3 result 
is simply owing to the limited sample size. 

The correlation length is smaller in the duty cycle sce- 
nario than in the escape fraction scenario with ro ~ 
3/i Mpc, which is reasonable given the lower mean 
stellar mass of LAEs in this scenario. The sparse sam- 
pling also prohibits us from obtaining the CF signal at 
r < 1.5 ft. -1 Mpc in the duty cycle scenario. In order to 
measure the CF for the duty cycle scenario more reliably, 
we need a simulation box size of Lbox ^ 100/i _1 Mpc 



with a resolution comparable to that of the D5 run. 
This should become possible in the near future thanks 
to rapidly increasing computing resources. 

11.1. Bias of LAEs 
The lower panels of Fig. [5] show the bias relative to the 
clustering of the mass, which is defined as b = \J £ ga i/£dm- 
We compute the correlation function of dark matter by 
randomly sampling 200,000 dark matter particles in the 
G6 run. Using the D5 run yields very similar results on 
scales of r = 1 — 8/i _1 Mpc. In both scenarios, the bias 
is a slowly decreasing function with increasing distance. 
Even though the value of ro is greater at z = 3, the 
bias relative to the dark matter is greater at z — 6, be- 
cause the growth in dark matter structure significantly 
increases the normalization of the dark matter CF from 
z = 6 to 3. 

In the escape fraction scenario (Fig. [5^), we find b ~ 3.5 
(5.0) at r = 1.5 — 10 h~ 1 Mpc for z = 3 (6), crossing the 
above value at r = 4 — 5/i _1 Mpc. At smaller scales 
(r < 1.5/i _1 Mpc), the bias increases up to b ~ 6 (9) 
at z = 3 (6). This could owe to the exce ss clustering 
of gala xies on small scales as discovered by lOuchi et al.l 
(|2005a| ). although Ouchi's data at z = 4.0 suggest that 
this increase in bias occurs at r < 0.2 ft, -1 Mpc. The 
increase of clustering on small-scales can be ascribed to 
the substructures within each halo ("one- halo" term). 
It is possible that the simulation is still lacking some 
physics or resolution to capture the correct scale for this 
transition. On large scales, the correlation function of 
dark matter seems to turn down, and at z = 6 the bias 
somewhat increases atr > 15/i _1 Mpc, which is probably 
a box-size effect. 
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In the duty cycle scenario (Fig. [8p), we use both the 
power-law fits and the actual simulation results to com- 
pute the bias. The calculation using the power-law yields 
monotonically declining functions with increasing dis- 
tance. At z = 3 the bias decreases from 4.5 to 1.5, and at 
z = 6 from 5.3 to 3.7 with increasing distance. The cal- 
culation using the direct simulation CF shows more noisy 
behavior, wiggling around the power-law result. At z = 3 
the bias decreases from 3.8 to 1.4 at r = 1.6 — 8 h^ 1 Mpc. 
At z — 6, the wiggle is smaller and the long-short-dashed 
line agrees well with the power-law result, yielding b ~ 4 
at r — 1.5 — 7 h^ 1 Mpc, which is contrasted with b ~ 5 
in the escape fraction scenario. 

The comparison to some of the observational esti- 
mates yields somewhat mixed res ults, but in g enera l 
support the duty cycle scenario. IKovac et all (j2007l ) 
found ro = 4.61 ± 0.6 ft. -1 Mpc (taking the contami- 
nation by randomly distributed objects into account) 
and b ~ 3.7 for the LAE s at z ~ 4.5 i n the LALA 
survey (jRhoads et al.ll2000D . lOuchi et al.l (120031) found 
r = 3 .5 ± 0.3 h' 1 Mpc for z = 4.86 LAEs. IKovac et all 
(|2007l ) pointed out that Ouchi's maximum permitted 
value would be r = 4.5 ± 0.4 h' 1 Mpc when the 20% 
contamination by low-z galaxies (Shimas aku et al.ll2004h 
is taken into account. Our results for the escape fraction 
scenario at z — 3 & 6 nicely bracket Kovac et al.'s results 
at z ~ 4.5. The bias values in the duty cycle scenario 
brackets the Kovac et al.'s result, but the value of r is 
lower than theirs or at the lower edge of Ouchi et al.'s 
e stimate. 

iGawiser et ail (|2007t) reported r = 3.6^ Mpc and 

b = l.TtoA for LAEs at z = 3A - 0ur CF results at z = 3 
for the duty cycle scenario are in good agreement with 
Gawiser et al.'s estimates. This is consistent with the 
nice agreement between our simulations and the observed 
data seen in Fig. [5] for the duty cycle scenario at z = 3. 
Our results from the escape fraction scenario at z — 3 do 
not agree with Gawiser et al.'s data. At z = 6, our result 
of b = 3.7 at r ~ 8ft. -1 Mpc in the duty c ycle scenario 
agrees well with b — 3.4 ± 1.8 derived by I Ouchi et al.l 
(|2005bf) for LAEs at z ~ 5.7. These comparisons again 
support the duty cycle scenario. 

12. COSMIC VARIANCE 

It is clearly difficult to conduct a deep, wide-field sur- 
vey of high-z galaxies and obtain a statistically represen- 
tative sample of galaxies in a large volume of space. If 
the survey area is too small, then the observed sample 
may not be representative of the total population owing 
to the large-scale structure of the Universe, and the esti- 
mated LF would scatter around the true LF. This is one 
manifestation of the so-called "cosmic variance" . 

In order to estimate the cosmic variance in our results, 
we use eight sub- volumes of (45 x 45 x 44/i _1 Mpc) 3 in 
the G6 run and derive the LF from each sub- volume (Fig- 
ure^. For our adopted flat A cosmology, a viewing angle 
of 1 degree corresponds to about comoving 100 /i -1 Mpc 
at z — 6, so the above subvolume corresponds to a field 
with FoV = (0.45 deg) 2 rj 0.2 deg 2 . The thickness of 
44/T" 1 Mpc was chosen to match the data of lOuchi et al.l 
(|2007l Fig. 18), but the exact value is not so important, 
as we would obtain a similar result even if we adopted 
(50 h^ 1 Mpc) 3 subvolumes. Here we used fi, ya = 0.15 as 
we did in Fig. |4j 



Z = 6 
FoV = 0.2 deg 2 




Fig. 9. — Lyct luminosity function at z = 6 measured in the eight 
subvolumes of (45 X 45 X 44 h~ 1 Mpc) 3 in the G6 run, corresponding 
to a 0.2 deg 2 field of view. The error bars show Poisson errors in 
each subvolume, slightly offset from each other to avoid ove rlap. 
The red filled triangle data points are from lOuchi et alj(|2007h . and 
the yellow shade shows the variance of their five 0.2 deg 2 fields. 

The large scatter of LFs seen in Fig. [9] owes to cosmic 
variance, and it clearly exceeds Poisson errors shown by 
the error bars. T he red filled triangle is the data by 
I Ouchi etTa l. (2007) and the yellow shade is the variance 
of their data from five 0.2 deg fields. The cosmic vari- 
ance we find is consistent with t he field-to-field variation 
observe d by lOuchi et al. ( 2007 . Fig. 18), as well as the 
data of IShimasaku et all ( 2004 ) that shows very differ- 
ent distribution of LAEs at z = 4.79 and 4. 86, separated 
by ~ 40 h~ Y Mpc in the same FoV (see also lHu fc Cowid 
2006). At log^Lya = 43.0, the cosmic variance in the 
vertical direction is almost an order of magnitude in 
dn/dlogLLy a , and ~ 0.6dex in the horizontal direction 
(i.e., logLLya)- At logLLya > 43.5, the Poisson error bar 
is large, because one object in the above subvolume corre- 
sponds to the data point of dn/dlogLL yQ = — 4.72+^ 30 . 
Our result shows that a survey field of > 1 deg 2 is nec- 
essary to obtain a reliable estimate of Lya & UV LFs at 
z = 6. 

13. CONCLUSIONS & DISCUSSIONS 

Encouraged by the earlier successes in describing the 
properties of LBGs at z = 3 — 6 usin g cosmologic al SPH 
simulations (jNagamine et al.l l2004ch iNight et al1l2006h . 
we considered two simple scenarios to explain the lumi- 
nosity functions of LAEs at z — 3 and 6. These scenar- 
ios are very simple, but should capture the two extreme 
cases of LAE properties. The true physical nature of 
LAEs could be a combination of the two scenarios. We 
summarize the parameters associated with the two sce- 
narios in Table [2] The main conclusions of our work can 
be summarized as follows: 

• In our simulations, star formation becomes progres- 
sively less efficient from z = 6 to 3, especially for 
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the massive galaxies with M* > 10 10 M Q . This 
evolution in the SF efficiency largely compensates 
the growth in the galaxy mass function, yielding 
little evolution in the rest-frame UV LF. We pro- 
vide a fitting formula for the SFR-M* relationship 
in Eq. (J). The SFR function (Fig. 1} does not 
evolve very much from z — 6 to 3. 

We investigate two simple scenarios to character- 
ize the properties of LAEs: the "escape fraction!' 
scenario and the "duty cycle" scenario. By match- 
ing the simulated Lya LF of LAEs to the observed 
one bv lOuchi et al.( (|2007f) . we find that the effec- 
tive escape fraction of Lya photons is fhya = 0.1 
(0.15) for z = 3 (6), where the effect of IGM at- 
tenuation is included. In the case of the duty cycle 
scenario, we find that the duty cycle parameter is 
Cduty = 0.07 (0.2) for z = 3 (6), after correcting 
the observed LF for the IGM attenuation effect. If 
we do not correct for IGM attenuation, we obtain 
Cduty = 0.06 for both z = 3 and 6. 

In both scenarios, we find that /dust ~ 0.1 at z = 3 
and /dust ~ 0.2 — 0.4 at z = 6. In detail, if we 
correct for the IGM attenuation effect with /igm = 
0.82 (0.52) and assume f™ = 0.06 (0.2) at z = 3 
(6), then we obtain /dust = 0.13 (0.36) at z = 3 (6), 
respectively. If we instead adopt /igm = 1.0 and 
the same values of /™\ then we find /dust =0.11 
(0.19) for z = 3 (6). 

The mean stellar mass of LAEs is lower in the duty 
cycle scenario than in the escape fraction scenario, 
as summarized in Tabled This implies lower mean 
values of BH mass, SFR and metallicity, as well as 
a lower amplitude of the correlation function in the 
duty cycle scenario. The mean stellar mass of LAEs 
is higher at z = 3 than at z = 6, as expected in the 
hierarchical structure formation model. 

In the escape fraction scenario, the auto-correlation 
function of LAEs is similar to that of LBGs with 
ro = 4 — 6/i _1 Mpc and 7 ~ 1.7 for both z = 3 
and 6. The corresponding galaxy-dark matter bias 
is b ~ 3.5 (5.0) at r = 1.5 - lO/i^Mpc for z = 3 
(6). The se results bracket t he observational esti- 
mates by iKovac et al.1 (|2007| ) at z ~ 4.5. On the 
other hand, the CF in the duty cycle scenario has 
a shorter correlation length of ro ~ 3 hr 1 Mpc and 
the bias of b = 1.6 - 4.6 [z = 3) & 4 {z = 6), in 
go od agreement with t he observational estimates 
by iGawiser et al.l (|2007l ) at z ~ 3. In both scenar- 
ios, the bias parameter decreases with increasing 
distance. 

We find that the effect of cosmic variance on LF 
estimates can be quite strong. The Lya LFs at 
z = 6 measured in 8 different fields of ~ 0.2 deg 2 
have a variance of ~ 1.0 dex in dn/dlogLLya at 
log-LLy Q = 43.0, and ~ 0.6 dex in the horizontal 
direction (i.e., logLLya)- This result suggests that 
a survey field of > (1 deg) 2 is necessary to measure 
the Lya and UV LFs at z = 6 reliably. 



• We find that the duty cycle scenario is preferred 
over the escape fraction scenario through the com- 
parisons with various observational data, includ- 
ing M*, SFR, Lya & UV LFs, correlation func- 
tions, and bias relative to the dark matter dis- 
tribution. In particular, the duty cycle scenario 
succeeds in explaining the the rest-frame UV LFs 
of both LBGs and LAEs at z = 3 & 6, provided 
E(B — V) ~ 0.15 for both population. How- 
ever, if the extinction of UV photons in LAEs 
is systematically lower than that of LBGs with 
E(B — V) < 0.05 as the recent observatio ns suggest 
(|Ouchi et alJl2007L iGronwall et al.ll2007f l. the n our 
simple models (with uniform f^ ya or Cduty) fail to 
explain the UV LFs of LAEs and LBGs simulta- 
neously. In this case, the escape fraction scenario 
completely breaks down, because it would predict 
a brighter UV LF of LAEs at z = 3, contrary to 
the observation. 

In the case of duty cycle scenario, additional model- 
ing would be necessary to accommodate lower extinction 
for LAEs, such as variable extinction and treating LBGs 
and LAEs as separate populations. The assumption of 
linear SFR - L^ ya relationship (Eq. [2|) and uniform /dust 
& E(B — V) imply that galaxies with high SFR would 
appear as bright LBGs with high L-^ ya . If we further as- 
sume th at E(B—V) is positive ly correlated with metallic- 
ity as in lFinlator et al.l (|2006h . then the mass- metallicity 
relationship implies that the bright LBGs would be rela- 
tively massive galaxies with high extinction. Combined, 
this means that high L^ ya sources have high extinction, 
which is the opposite of what is found observationally 
( Shap lev et all 12003 lOuchi et al.l I2007L IGronwall et all 
l2007t IGawiser et HTf l2007l )~ Therefore, to assign lower 
extinction for LAEs, one must invoke some other physi- 
cal processes, such as galactic outflows, that could tem- 
porarily reverse the above relationship between L\^ ya and 
E{B-V). 

If the duty cycle scenario is really the correct one, then 
it could explain both Lya and UV LFs of LAEs as fol- 
lows: in the universe as early as z = 6, LBGs are ac- 
tively forming stars, and as much as 20% (Cduty = 0.2) 
of them would appear as LAEs (right panel of Fig. [7]). 
If our simulations overpredict the UV LF of LBGs at 
2 = 6, then this fraction could increase even further. By 
the time the universe has evolved to z = 3, star for- 
mation becomes less efficient owing to the lower mean 
density and less supply of inf ailing gas than at z = 6. 
The extinction of Lya photons by dust also becomes 
stronger at z = 3 (/dust ~ 0.1). Only a small frac- 
tion (Cduty = 0.07) of LBGs now appear as LAEs. 
Ou r values of Cdutv nicely br acket the result obtained 
bv iMalhotra fc Rhoadsl (|2002l C du t y = 0.075 - 0.15 at 
z = 4.5). If the star formation time-scale for the z — 3 
LBGs is ~ 200 Myr, then the life-time of LAEs at z = 3 
would be ~ 14 Myr according to the above duty cycle 
parameter. Ages between several Myr to 100 Myr are 
supported by re ce nt observations (IGawiser et alj 120061 : 
lLai et al.|[2007aHbl : iPirzkal et al J 120071 ). We remark that 
sporadic star formation histories are automatically in- 
cluded in our results, sin ce they occur generical ly in our 
dynamic simulations (see iNagamine et al.ll2005T ). There- 
fore, the values of Cduty in this paper characterize the 
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duty cycle on top of the intermittency of the SF activity. the two work are consistent with each other. 



Because we adopted a single relationship between SFR 
and Lya luminosity (Eq. [2]) , we are essentially assuming 
that the equivalent width (EW) distribution of the Lya 
emission line is a ^-function at ~ 70 A. 6 This is obviously 
oversimplified, given that the actually obser ved data has 
a fair ly wide distribution (e.g., Fig. 23 of lOuchi et al.1 
12001 . One can easily accommodate a model that as- 
sumes a certain distribution of EW by, e.g., linking it to 
the scatter in the dust and metallicity distribution , and 
observe its effect on the LF. iKobavashi et al.l (|2007h per- 
formed such an exercise using a semi-analytic model of 
galaxy formation, and showed that the bright-end of the 
Lya LF can be described better with such a treatment. 
Currently it is not clear whether our simulation is over- 
predicting the Lya LF at the bright-end, because there 
are no reliable data at logLLya > 43.5 (see Fig. |4|). If 
future observations show that there is a sharp exponen- 
tial cutoff at logLLya > 43.5, we would be able to make 
stronger arguments on the possible connections between 
feedback effects and the bright-end of the Lya LF. We 
note that our current simulations already include the ef- 
fect of SN feedback and galactic outflows as described in 
§ [2] Given these uncertainties, we decided not to invoke 
an additional model for the EW distribution on top of 
what has been modeled in the simulation. 



IDave et al.l (|2006f) computed the Lya LF using cosmo- 
logical SPH simulation in a similar fashion to our present 
work. They assume d Lya emission of 2.44 x 10 42 ergs -1 
per SFR [M yr" 1 ] (|Schaereiil2003l) and took the metal- 
licity variation into account. They a rgued that they ob - 
tain good agreement with the data of lSantos et al.1 (|2004D 
at logiLya = [40.5,42.5] if th ey assume fLya = -02. 
However, the LF data points of (Santos et~aU(|200l are 
based on only a few objects and their normalization 
seems to be too low t o smo othly connect with the data 
points of lOuchi et al.1 (|2007f ) at logLLya > 42, as shown 
by the yellow shade in Fig. 0J Our simulations suggest 
that the LF at Li, ya = [40 5,42.5] woul d be h igher by a 
factor of ~ 5 than that of lSantos et al.l (l2004h whe n we 
match our results to the data o f lOuchi et al.l (120071) . We 
also point out that, according to wh at we discu s sed in 
§ [121 the simulation box size used by IDave et"aLl (|2006l ) 
(comoving 33/i _1 Mpc) was too small to obtain a reliable 
LF at the bright-end. Had they used a simulation with 
a l arger box size and matched their results to the data 
by lOuchi et alj (|2007f ). we expect that they would have 
obtained a value of fhya similar to ours. Their galaxy 
correlation function at z = 6 is also lower than ours, 
probably owing to the same reason of a small box size. 
We note, however, that the galaxy stellar mass function 
(in the mass range that they simulated) and the bias in 

6 We calculated the EW as follows: EW = L Lya /L x (1216A) = 
^Lya/([c/A 2 ]L 1/ [1216A]) = 70 A, where we used Lt vc , 

1 erg s" 1 per SFR from Eq. J2J|, and L„(1216A) = L„(1500A) = 
7xl0 27 ergs" 1 Hz -1 per SFR for a Salpeter IMF with mass limits 



The value of .fhv a = 0-02 assumed by both 
IDave et al.l (|2006ft and iLe Delliou et~aTl (|2006D is signif- 
icantly smaller than our value of fhya — 0.1. It is not 
easy to understand the source of this difference unless 
we compare the SFR function (Fig. [2]) and the effect of 
dust a ttenuation in each model in detail. At least in the 
case of lDave et all (|2006l ). they may have underestimated 
the value of fhya by adjusting the results from a sma ll 
simulation box size to the data oflSantos et ail (|2004l) . 
On the other hand, IKobavashi et al.l ~( 20071 ) argued for 
a much larger value for the escape fraction, / csc ~ 0.8, 
by incorporating the effects of outflows from starburst 
galaxies. As we described earlier, our simple escape frac- 
tion model fails to reproduce the UV LF of LAEs if their 
extinction is systematically lower than that of LBGs. It 
would be interesting to see what the above semi-analytic 
models predict on the joint constraint of the UV LF of 
LAEs and LBGs at z = 3 - 6. 

Another source of uncertainty in the present work is 
the IMF. It is possible that the IMF changes as a f unction 
of re dshift, metallicity, and environment (e.g., iLarsonl 
Il998| ) . Recently there have been several suggestions that 
some of the observational data can be better explained if 
the IMF be comes "top-heavy" or "bottom-light" towards 
high- z fe.g. jBaugh et al.ll2005tlChary||2007tlvan Dokkuml 
12007ft . It would also mitigate the apparent discrepancy 
between the data on extinction corrected SFR density 
and stellar mass densi t y (e.g.. | Nagamine et al.l 2004a; 
iNaeamine et all l2006al : iFardal et all 120071 : iDavel 12007 : 
IWilkins et al.ll2008l ). If indeed the IMF is top-heavy at 
high-z, then the convers ion factor in E q. @ could in- 
crease by a factor of ~ 2 (Schacrcr 2003|), and our intrin- 
sic Lya LF would be brighter by the same factor. In this 
case, the values of fhya and Cduty should be decreased 
by similar factors. The top-heavy IMF would also dis- 
rupt the agreement of the UV LF of LBGs at z = 3 
between simulations and the observed data, and exac- 
erbate the discrepancy between the simulation and the 
observed data at z = 6. As far as the Lya & UV LFs 
of LAEs are concerned, we do not see any strong reasons 
to invoke a top-heavy IMF based on our comparisons, 
although we cannot rule out the possibility of top-heavy 
IMF given various uncertainties discussed in this paper. 
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TABLE 2 
Summary of Two Scenarios 



Parameter Escape Fraction Scenario Duty Cycle Scenario 





2 = 3 


2 = 6 


2 = 3 


2 = 6 


/tya 


0.10 


0.15 


— 


— 


/iGM b 


0.82 (1.0) 


0.52 (1.0) 


0.82 (1.0) 


0.52 (1.0) 


j"ion c 
J esc 


0.06 


0.20 






/dust d 
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0.36 (0.19) 






c USt e 

duty 






0.07 (0.06) 


0.2 (0.06) 


log(Af*/Af Q ) f 


[10.2, 12.1] 


[9.3, 11.0] 


[9.2,11.3] 


[8.6,10.3] 


(M*/M @ ) a 


2.5 x 10 10 


1.9 x 10 9 
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p-k fraction h 
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0.11 


0.42 


0.29 


log(M BH /M Q ) i 
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[6.9,8.7] 


[6.7,8.8] 
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(M BH /M > J 
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(Z/Z ) m 
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0.17 


0.21 


0.11 


ro n 


5.5 


4.4 


(3.2) 


3.1 


7 P 


1.67 


1.68 


(2.30) 


1.49 


& q 


3.5 (< 6.0) 


5.0 (< 8.9) 


(1.6-4.6) 


4.0 



Note. — Summary of various parameter values in the two scenarios discussed in 
this paper. 

a: Effective escape fraction of Lyo photons including the effect of IGM attenuation 

(§E3J. 

b: IGM attenuation factor computed by following Madau (1995) prescription (§ 15.11 1, 
The case for /igm = 1-0 is given in parenthesis. 

c: Escape fraction of ionizing photons taken from llnoue et al.l f200^ . Fig. 3) (§1511. 
d: Extinction of Lya photons by the local dust (§ 15.111 . The case for /igm = 1-0 is 
given in parenthesis. 

e: Duty cycl e (or fractional life time) of LAEs after correcting for the IGM attenuation 
effect (§ I5.2J| , The case for /igm = 1-0 is given in parenthesis. 
/: Stellar mass range of LAEs with logLLya = 42 — 44 (§|6j. 
g: Mean stellar mass of LAEs with logLL ya = 42 — 44 (§(6j. 

h: Fraction of stellar mass density contained in LAEs with logLLya =42 — 44 (§0. 
i: Mass range of black holes hosted by LAEs with logLLya = 42 — 44 (§0. 
j: Mean mass of BHs hosted by LAEs with logLLya 

= 42 - 44 (§12). 

k: Stellar metallicity range LAEs with logI/ L y Q =42 — 44 (§0. 
m: Mean stellar metallicity of LAEs with logLLya = 42 — 44 (§0. 
n: Correlation length of LAEs (§ 1111 1. 
p: Slope of the auto-correlation function of LAEs (§ I1H . 

q: Bias parameter of LAEs against dark matter distribution. For the escape fraction 
scenario, the maximum values at r ~ 0.9 /i — 1 Mpc is given in the parenthesis. The 
result for 2 = 3 duty cycle scenario is given in parenthesis, because they are somewhat 
unreliable owing to small sample size. 



Baugh, C. M., Lacey, C. G., Frenk, C. S., Granato, G. L., Silva, 
L., Bressan, A., Benson, A. J., & Cole, S. 2005, MNRAS, 356, 
1191 

Beckwith, S. V. W., Stiavelli, M., Koekemoer, A. M., Caldwell, 

J. A. R., Ferguson, H. C, Hook, R., Lucas, R. A., Bergeron, 

L. E., et al. 2006, AJ, 132, 1729 
Bouwens, R. J., Illingworth, G. D., Blakeslee, J. P., & Franx, M. 

2006, ApJ, 653, 53 
Bouwens, R. J., Illingworth, G. D., Franx, M., & Ford, H. 2007, 

ArXiv e-prints, 707 
Bouwens, R. J., Illingworth, G. D., Thompson, R. I., Blakeslee, 

J. P., Dickinson, M. E., Broadhurst, T. J., Eisenstein, D. J., 

Fan, X., et al. 2004, ApJL, 606, L25 
Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 
Bunker, A. J., Stanway, E. R., Ellis, R. S., & McMahon, R. G. 

2004, MNRAS, 355, 374 
Chary, R.-R. 2007, ArXiv e-prints, 712 

Chen, H.-W., Prochaska, J. X., & Gnedin, N. Y. 2007, ApJ, 667, 
L125 

Dave, R. 2007, ArXiv e-prints, 710 

Dave, R., Finlator, K., & Oppenheimer, B. D. 2006, MNRAS, 370, 
273 

Dave, R., Hernquist, L., Katz, N., & Weinberg, D. H. 1999, ApJ, 
511, 521 



Dawson, S., Rhoads, J. E., Malhotra, S., Stern, D., Dey, A., 
Spinrad, H., Jannuzi, B. T., Wang, J., & Landes, E. 2004, ApJ, 
617, 707 

Di Matteo, T., Colberg, J., Springel, V., Hernquist, L., & Sijacki, 

D. 2007, ArXiv e-prints (0705.2269), 705 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 
604 

Dickinson, M., Stern, D., Giavalisco, M., Ferguson, H. C, 
Tsvetanov, Z., Chornock, R., Cristiani, S., Dawson, S., Dey, A., 
Filippenko, A. V., Moustakas, L. A., Nonino, M., Papovich, C, 
Ravindranath, S., Riess, A., Rosati, P., Spinrad, H., & Vanzclla, 

E. 2004, ApJ, 600, L99 

Dow-Hygelund, C. C, Holden, B. P., Bouwens, R. J., Illingworth, 
G. D., van der Wei, A., Franx, M., van Dokkum, P. G., Ford, H., 
et al. 2007, ApJ, 660, 47 

Drory, N., Salvato, M., Gabasch, A., Bender, R., Hopp, U., Fculncr, 

C, & Pannella, M. 2005, ApJ, 619, L131 

Fardal, M. A., Katz, N., Gardner, J. P., Hernquist, L., Weinberg, 

D. H., & Dave, R. 2001, ApJ, 562, 605 

Fardal, M. A., Katz, N., Weinberg, D. H., & Dave, R. 2007, 

MNRAS, 379, 985 
Finlator, K., Dave, R., Papovich, C, k, Hernquist, L. 2006, ApJ, 

639, 672 



16 



Nagamine et al. 



Fontana, A., Salimbeni, S., Grazian, A., Giallongo, E., Pentcricci, 
L., Nonino, M., Fontanot, F., Mcnci, N., Monaco, P., Cristiani, 
S., Vanzella, E., de Santis, C, & Gallozzi, S. 2006, A&A, 459, 
745 

Furlanctto, S. R., Schaye, J., Springel, V., & Hernquist, L. 2005, 
ApJ, 622, 7 

Gawiscr, E., Francke, H., Lai, K., Schawinski, K., Gronwall, C, 
Ciardullo, R., Quadri, R., Orsi, A., et al. 2007, ArXiv c-prints, 
710 

Gawiscr, E., van Dokkum, P. G., Gronwall, C, Ciardullo, R., Blanc, 

G. A., Castandcr, F. J., Feldmeier, J., Francke, H., et al. 2006, 
ApJ, 642, L13 

Gronwall, C, Ciardullo, R., Hickey, T., Gawiser, E., Feldmeier, 

J. J., van Dokkum, P. G., Urry, C. M., Herrera, D., et al. 2007, 

ApJ, 667, 79 
Haiman, Z. & Hui, L. 2001, ApJ, 547, 27 
Haiman, Z., Spaans, M., & Quataert, E. 2000, ApJ, 537, L5 
Hcckman, T. M. 2001, in Astronomical Society of the Pacific 

Conference Series, Vol. 240, Gas and Galaxy Evolution, ed. J. E. 

Hibbard, M. Rupen, & J. H. van Gorkom, 345 
Hopkins, P. F., Cox, T. J., Kcres, D., & Hernquist, L. 2007a, ArXiv 

e-prints, 706 

Hopkins, P. F., Hernquist, L., Cox, T. J., & Kcres, D. 2007b, ArXiv 
e-prints, 706 

Hopkins, P. F., Hernquist, L., Cox, T. J., Robertson, B., & Krausc, 

E. 2007c, ApJ, 669, 45 
Hu, E. M. & Cowie, L. L. 2006, Nature, 440, 1145 
Hu, E. M., Cowie, L. L., Capak, P., McMahon, R. G., Hayashino, 

T., & Komiyama, Y. 2004, AJ, 127, 563 
Hu, E. M. & McMahon, R. G. 1996, Nature, 382, 231 
Inoue, A. K., Iwata, I., & Deharveng, J.-M. 2006, MNRAS, 371, 

LI 

Iwata, I., Ohta, K., Tamura, N., Akiyama, M., Aoki, K., Ando, M., 

Kiuchi, G., & Sawicki, M. 2007, MNRAS, 376, 1557 
Kashikawa, N., Kitayama, T., Doi, M., Misawa, T., Komiyama, Y., 

& Ota, K. 2007, ArXiv e-prints, 704 
Kashikawa, N., Shimasaku, K., Malkan, M. A., Doi, M., Matsuda, 

Y., Ouchi, M., Taniguchi, Y., Ly, C, et al. 2006, ApJ, 648, 7 
Katz, N., Weinberg, D. H., & Hernquist, L. 1996, ApJS, 105, 19 
Keel, W. C, Cohen, S. H., Windhorst, R. A., & Waddington, I. 

1999, AJ, 118, 2547 
Kennicutt, R. C. J. 1998, ApJ, 498, 541 

Kobayashi, M. A. R., Totani, T., & Nagashima, M. 2007, ArXiv 
e-prints, 705 

Kormendy, J. & Gcbhardt, K. 2001, in AIP Conf. Proa, 20th Texas 
Symposium on Rclativistic Astrophysics, ed. J. C. Wheeler & 

H. Martel, Vol. 586 (New York: Am. Inst. Phys.), 363 

Kovac, K., Somerville, R. S., Rhoads, J. E., Malhotra, S., & Wang, 

J. 2007, ArXiv e-prints, 706 
Lai, K., Huang, J.-S., Fazio, G., Cowie, L. L., Hu, E. M., & Kakazu, 

Y. 2007a, ApJ, 655, 704 
Lai, K., Huang, J.-S., Fazio, G., Gawiser, E., Ciardullo, R., Damen, 

M., Franx, M., Gronwall, C, Labbe, I., Magdis, G., & van 

Dokkum, P. 2007b, ArXiv e-prints, 710 
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64 
Larson, R. B. 1998, MNRAS, 301, 569 

Le Delliou, M., Lacey, C, Baugh, C. M., Guiderdoni, B., Bacon, 
R., Courtois, H., Sousbie, T., & Morris, S. L. 2005, MNRAS, 
357, Lll 

Le Delliou, M., Laccy, C. G., Baugh, C. M., & Morris, S. L. 2006, 

MNRAS, 365, 712 
Leitherer, C, Schaerer, D., Goldader, J. D., Dclgado, R. M. G., 

Robert, C, Kunc, D. F., dc Mcllo, D. F., Devost, D., & Heckman, 

T. M. 1999, ApJS, 123, 3 
Li, Y., Hernquist, L., Robertson, B., Cox, T. J., Hopkins, P. F., 

Springel, V., Gao, L., Di Matteo, T., Zentner, A. R., Jenkins, 

A., & Yoshida, N. 2007, ApJ, 665, 187 
Madau, P. 1995, ApJ, 441, 18 

Madau, P., Pozzetti, L., & Dickinson, M. 1998, ApJ, 498, 106 
Malhotra, S. & Rhoads, J. E. 2002, ApJ, 565, L71 
— . 2004, ApJ, 617, L5 

Malhotra, S., Rhoads, J. E., Pirzkal, N., Haiman, Z., Xu, C, Daddi, 
E., Yan, H., Bergeron, L. E., et al. 2005, ApJ, 626, 666 

Matsuda, Y., Yamada, T., Hayashino, T., Tamura, H., Yamauchi, 
R., Ajiki, M., Fujita, S. S., Murayama, T., et al. 2004, AJ, 128, 
569 

McQuinn, M., Hernquist, L., Zaldarriaga, M., & Dutta, S. 2007, 
MNRAS, 381, 75 



Mori, M. & Umemura, M. 2006, Nature, 440, 644 

Murayama, T., Taniguchi, Y., Scoville, N. Z., Ajiki, M., Sanders, 
D. B., Mobashcr, B., Aussel, H., Capak, P., Kockcmoer, A., 
Shioya, Y., Nagao, T., Carilli, C, Ellis, R. S., Garilli, B., 
Giavalisco, M., Kitzbichlcr, M. G., Le Fevre, O., Maccagni, D., 
Schinnerer, E., Smolcic, V., Tribiano, S., Cimatti, A., Komiyama, 
Y., Miyazaki, S., Sasaki, S. S., Koda, J., & Karoji, H. 2007, ApJS, 
172, 523 

Nagamine, K., Cen, R., Hernquist, L., Ostriker, J. P., & Springel, 

V. 2004a, ApJ, 610, 45 
— . 2005, ApJ, 627, 608 

Nagamine, K., Ostriker, J. P., Fukugita, M., & Cen, R. 2006a, ApJ, 
653, 881 

Nagamine, K., Springel, V., & Hernquist, L. 2004b, MNRAS, 348, 
435 

Nagamine, K., Springel, V., Hernquist, L., & Machacck, M. 2004c, 

MNRAS, 350, 385 
Nagamine, K., Wolfe, A. M., & Hernquist, L. 2006b, ApJ, 647, 60 
Nagashima, M., Yahagi, H., Enoki, M., Yoshii, Y., & Gouda, N. 

2005, ApJ, 634, 26 
Night, C, Nagamine, K., Springel, V., & Hernquist, L. 2006, 

MNRAS, 366, 705 
Oesch, P. A., Stiavelli, M., Carollo, C. M., Bergeron, L. E., 

Koekemoer, A. M., Lucas, R. A., Pavlovsky, C. M., ct al. 2007, 

ArXiv e-prints, 706 
Ouchi, M., Hamana, T., Shimasaku, K., Yamada, T., Akiyama, M., 

Kashikawa, N., Yoshida, M., Aoki, K., Iye, M., Saito, T., Sasaki, 

T., Simpson, C, & Yoshida, M. 2005a, ApJ, 635, L117 
Ouchi, M., Shimasaku, K., Akiyama, M., Sekiguchi, K., Furusawa, 

H., Okamura, S., Kashikawa, N., Iye, M., Kodama, T., Saito, T., 

Sasaki, T., Simpson, C, Takata, T., Yamada, T., Yamanoi, H., 

Yoshida, M., & Yoshida, M. 2005b, ApJ, 620, LI 
Ouchi, M., Shimasaku, K., Akiyama, M., Simpson, C, Saito, T., 

Ueda, Y., Furusawa, H., Sekiguchi, K., Yamada, T., Kodama, 

T., Kashikawa, N., Okamura, S., Iye, M., Takata, T., Yoshida, 

M., & Yoshida, M. 2007, ArXiv e-prints, 707 
Ouchi, M., Shimasaku, K., Furusawa, H., Miyazaki, M., Doi, 

M., Hamabc, M., Hayashino, T., Kimura, M., Kodaira, K., 

Komiyama, Y., et al. 2003a, ApJ, 582, 60 
Ouchi, M., Shimasaku, K., Furusawa, H., Miyazaki, M., Doi, M., 

Hamabe, M., Hayashino, T., Kimura, M., et al. 2003, ApJ, 582, 

60 

Ouchi, M., Shimasaku, K., Okamura, S., Furusawa, H., Kashikawa, 
N., Ota, K., Doi, M., Hamabe, M., Kimura, M., Komiyama, Y., 
Miyazaki, M., Miyazaki, S., Nakata, F., Sekiguchi, M., Yagi, M., 
& Yasuda, N. 2004a, ApJ, 611, 660 

— . 2004b, ApJ, 611, 685 

Partridge, R. B. & Peebles, P. J. E. 1967, ApJ, 147, 868 

Peng, C. Y., Impey, C. D., Rix, H.-W., Kochanck, C. S., Keeton, 

C. R., Falco, E. E., Lehar, J., & McLeod, B. A. 2006, ApJ, 649, 

616 

Pentericci, L., Grazian, A., Fontana, A., Salimbeni, S., Santini, 

P., De Santis, C, Gallozzi, S., & Giallongo, E. 2007, ArXiv 

Astrophysics c-prints 
Perez-Gonzalez, P. G., Riekc, G. H., Villar, V., Barro, G., Blaylock, 

M., Egami, E., Gallego, J., Gil de Paz, A., Pascual, S., Zamorano, 

J., & Donley, J. L. 2007, ArXiv e-prints, 709 
Pcttini, M. 2004, in Cosmochemistry. The melting pot of the 

elements, ed. C. Esteban, R. Garcia Lopez, A. Hcrrcro, & F. a. -p. 

Sanchez, 257-298 
Pettini, M., Rix, S. A., Stcidcl, C. C, Adelberger, K. L., Hunt, 

M. P., & Shapley, A. E. 2002, ApJ, 569, 742 
Pettini, M., Shapley, A. E., Steidel, C. C, Cuby, J.-G., Dickinson, 

M., Moorwood, A. F. M., Adelberger, K. L., & Giavalisco, M. 

2001, ApJ, 554, 981 
Pirzkal, N, Malhotra, S., Rhoads, J. E., & Xu, C. 2007, ApJ, 667, 

49 

Rhoads, J. E., Malhotra, S., Dcy, A., Stern, D., Spinrad, H., & 

Jannuzi, B. T. 2000, ApJ, 545, L85 
Robertson, B., Hernquist, L., Cox, T. J., Di Matteo, T., Hopkins, 

P. F., Martini, P., & Springel, V. 2006, ApJ, 641, 90 
Saito, T., Shimasaku, K., Okamura, S., Ouchi, M., Akiyama, M., 

Yoshida, M., & Ueda, Y. 2007, ArXiv e-prints, 705 
Santos, M. R., Ellis, R. S., Kneib, J.-P., Richard, J., & Kuijken, K. 

2004, ApJ, 606, 683 
Schaerer, D. 2003, A&A, 397, 527 

Shapley, A. E., Stcidcl, C. C, Adelberger, K. L., Dickinson, M., 
Giavalisco, M., & Pettini, M. 2001, ApJ, 562, 95 



LAEs and LBGs at z = 3 - 6 



17 



Shapley, A. E., Steidel, C. C, Pettini, M., & Adclbergcr, K. L. 

2003, ApJ, 588, 65 
Shimasaku, K., Hayashino, T., Matsuda, Y., Ouchi, M., Ohta, K., 

Okamura, S., Tamura, H., Yamada, T., et al. 2004, ApJ, 605, 

L93 

Shimasaku, K., Kashikawa, N., Doi, M., Ly, C, Malkan, M. A., 
Matsuda, Y., Ouchi, M., Hayashino, T., et al. 2006, PASJ, 58, 
313 

Sijacki, D., Springel, V., Di Matteo, T., & Hernquist, L. 2007, 

ArXiv e-prints, 705 
Springel, V. 2005, MNRAS, 364, 1105 

Springel, V., Di Matteo, T., & Hernquist, L. 2005, ApJ, 620, L79 
Springel, V. & Hernquist, L. 2002, MNRAS, 333, 649 
— . 2003a, MNRAS, 339, 289 
— . 2003b, MNRAS, 339, 312 

Stanway, E. R., Bunker, A. J., McMahon, R. G., Ellis, R. S., Treu, 

T., & McCarthy, P. J. 2004, ApJ, 607, 704 
Steidel, C. C, Adelbcrgcr, K. L., Shapley, A. E., Pettini, M., 

Dickinson, M., & Giavalisco, M. 2000, ApJ, 532, 170 



Taniguchi, Y., Ajiki, M., Nagao, T., Shioya, Y., Murayama, T., 
Kashikawa, N., Kodaira, K., Kaifu, N., et al. 2005, PASJ, 57, 
165 

van Dokkum, P. 2007, ArXiv e-prints, 710 

Vanzclla, E., Cristiani, S., Dickinson, M., Kuntschncr, H., Nonino, 
M., Rettura, A., Rosati, P., Vernet, J., et al. 2006, A&A, 454, 
423 

Venemans, B. P., Rottgering, H. J. A., Miley, G. K., Kurk, J. D., 
De Breuck, C, Overzier, R. A., van Breugel, W. J. M., Carilli, 
C. L., et al. 2005, A&A, 431, 793 

Wilkins, S. M., Trentham, N., & Hopkins, A. M. 2008, ArXiv c- 
prints, 801 

Yan, H. & Windhorst, R. A. 2004, ApJ, 612, L93 
Yoshida, M., Shimasaku, K., Kashikawa, N., Ouchi, M., Okamura, 
S., Ajiki, M., Akiyama, M., Ando, H., et al. 2006, ApJ, 653, 988 



